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Abstract. Weakly nonlocal thermodynamic theories are critically revisited. 
A relocalized, irreversible thermodynamic theory of nonlocal phenomena is 
given, based on a modified form of the entropy current and new kind of internal 
variables, the so called current multipliers. The treatment is restricted to 
deal with nonlocality connected to dynamic thermodynamic variables. Several 
classical equations are derived, including Guycr-Krumhansl, Ginzburg-Landau 
and Cahn-Hilliard type equations. 



1. Introduction 

In the last decades there has been a continuous interest in developing generalized 
classical continuum theories to include memory and nonlocal effects. By memory 
effects one means non instantaneous interactions and actions in the continuum, that 
is the dynamics of the continuum is determined not only by the present state but 
also by the previous ones, the history of the continuum. By nonlocal effects one 
means that the interaction between representative volume elements is not local, 
the dynamics is determined not only by the closest neighbours, the influence of 
farther volume elements is not negligible. There are at least three basic methods to 
consider memory and nonlocal effects in the dynamic equations. One can construct 
hypercontinuum theories introducing new variables beside the space-time variables 
(see e.g. P] E])- In case of strongly nonlocal theories the constitutive functions 
are given by a space integral form directly calculating the effects of the farther 
neighborhood (see e.g. [HIE])- The most popular are the so called weakly nonlocal 
theories, where higher order space derivatives are introduced into the governing 
equations than is usual in classical approaches. The original governing equations, 
that one wants to generalize, can be the balance equations of classical continuum 
physics or some other equations, too. The choice of the method is more or less a 
matter of convenience and experience, but without doubt the last method, the con- 
struction of weakly nonlocal continuum theories is the most attractive and popular 
among them . 

The popularity of weakly nonlocal theories is partially based on the convenience 
and the familiar knowledge of differential equations. The other approaches are 
somehow too general to be easily understandable and manageable, and bear some 
basic controversies. First of all, in hypcrcontinua or in the strongly nonlocal ap- 
proach some requirements of a physical theory (e.g. objectivity, distinction between 
material and structural components) are not included and exploited. Moreover, the 
nonlocal effects (e.g. a kind of weakening of interactions) are to be connected to the 
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Second Law at least in a general heuristical sense. Seemingly there is no clear con- 
ceptual prescription in the different formulations of the Second Law that would give 
a firm method in hypcrcontinua or strongly nonlocal theories. On the other hand, 
weakly nonlocal theories have an approximative character. First we can consider 
derivatives with one order higher than in the traditional equations and we can check 
the effects experimentally and theoretically. More systematic considerations should 
give a series of weakly nonlocal equations as a kind of successive approximation. 

It is interesting and instructive to review briefly the structure of the theories 
describing memory effects and compare them with nonlocal theories. The general 
memory functional method of rational thermodynamics (strong history dependence) 
(see e.g. in j^j) is the counterpart of strongly nonlocal theories. The time related 
counterparts of weakly nonlocal theories (weak history dependence) are the models 
that use rate dependent constitutive relations in the basic balance laws. Both are 
well treated and used approaches in this respect. However, the last development 
(more than thirty years ago ) was the introduction of internal variables of different 
kind, and in this way one can understand that the balance structure and the Second 
Law can generate the rate dependence (see e.g. in jjj). This method preserves 
the clear physical structure of classical irreversible thermodynamics and gives a 
considerable flexibility in modelling physical phenomena. Let us recognize that 
using internal variables instead of time derivatives in the constitutive relations, the 
constitutive relations themselves are strictly instantaneous, they depend only on 
state variables of the material, therefore only on the present state. The internal 
variables are related to the instantaneous state of the material, they are 'local' (in 
time), therefore, with their help the memory effects are modelled with a local theory, 
we can say that they are relocalized. The next table summarizes the situation: 





Space 


Time 


Strongly nonlocal 


space integrals 


memory functionals 


Weakly nonlocal 


gradient dependent con- 
stitutive functions 


rate dependent constitu- 
tive relations 


Relocalized 


??? 


internal variables 



In this paper we investigate the possibility to fill the place of the question marks 
in the table, to relocalize weak nonlocalities by introducing new kind of internal 
variables in a way that is compatible with the Second Law and with traditional 
methods of irreversible thermodynamics. We can accomplish our task when the 
nonlocality is connected to the so called dynamic internal variables (dynamic de- 
grees of freedom), those that vanish in (local) thermodynamic equilibrium. 

The general strategy of weakly nonlocal approaches is to include space deriva- 
tives, gradients of field quantities in the different classical constitutive functions 
(therefore sometimes they are called 'gradient theories'). The justification of these 
generalizations is based on either checking the consequences and/or checking the 
requirements of basic physical principles as, for example, different forms of objectiv- 
ity or/and the Second Law. Several approaches are used in mechanics e.g. different 
gradient theories [3 EUTUHIIl C2 U3HHI , the virtual power considerations [T31HB] 
and the multifield theories ^j^j. Instead of detailed reviewing and criticising the 
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different weakly nonlocal approaches our general remark is, that they usually intro- 
duce new and disputable concepts that seem to be too special to serve as foundation 
of a general nonlocal thermodynamic theory. Moreover, they are not constructed 
as a general approach to deal with nonlocal extensions of the transport equations of 
classical irreversible thermodynamics, they usually give methods for special fields 
of continuum physics (heat conduction, thcrmomechanics, etc.). 

1.1. Requirements of the Second Law. Introducing gradient dependence in 
an ad-hoc way one can get equations with solutions contradicting with the Sec- 
ond Law. Even if we keep the equations we should exclude some initial conditions 
and solutions because of physical reasons (e.g. because they contradict the experi- 
ments) . We should and can avoid these difficulties applying a systematic approach, 
a constitutive theory. To formulate a constitutive theory for nonlocal effects we 
strongly need a clear formulation of the Second Law and a construction method to 
get equations that automatically fulfill this requirement. The fundamental method 
of continuum physics to handle the situation is to separate the material and the 
structural properties of the phenomena: one should use the basic balance equations 
of the considered material and introduce constitutive functions that correspond to 
the Second Law. The question arises: what kind of Second Law, which formulation 
is the most fruitful? In this paper we accept that there is an extensive entropy 
function s, that should have a non negative production for supply free processes, in 
pure dissipative systems. We will see that this conservative standpoint is sufficiently 
general for our purposes. We do not need to weaken the well founded theoretical 
basis of irreversible thermodynamics and to exploit a freedom of any more general 
approach, e.g. to introduce a nonextensive entropy function. 

In non-equilibrium thermodynamics there are three traditional methods to check 
the compatibility of constitutive function with the inequality of the Second Law. 

- Irreversible thermodynamics derives an expression of the entropy produc- 
tion and introduces thermodynamic forces and currents with linear or non- 
linear relations between them ^3 [50] . 

- Coleman-Noll procedure investigates more systematically the linear inde- 
pendence of the different possible forces and currents [22 ■ 

- Liu procedure looks for the form of the constitutive functions that ensures 
the validity of the entropy inequality with any solutions of the governing 
differential equations, without substituting them into the entropy balance 

E2 

The methods are ordered with growing exactness. The first method is construc- 
tive. One can find a more detailed treatment regarding the different formulations 
and exploitations of the Second Law inequality for example in |23l 1241 l2~5] . In the 
following considerations one of our basic aims is to combine the heuristic power of 
the traditional force-current systems with a mathematically clearer treatment. 

1.2. Discrete and continuum. Beyond the compatibility with the Second Law in 
all continuum thermodynamic theories is important and instructive to clarify the 
relation between discrete and continuum descriptions of physical systems. Some 
continuum presentations pronounce to be self standing continuum theories without 
any connection to thermodynamics of homogeneous, discrete systems. Homoge- 
neous bodies are considered as special cases of continua. However, the relation is 
far more difficult if we consider processes in homogeneous thermodynamic systems 
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as it was shown in |2f>| . Anyway, a conceptual classification gives the basic heuris- 
tics. The relation of continuum and discrete systems is based on continuization 
principles. Classical irreversible thermodynamics is based on the principle of local 
equilibrium |19l I27| . i.e. in that theory we suppose that 

- physical quantities defined in ordinary (equilibrium) thermodynamics of 
discrete systems and the relations between them are valid in a continuum 
theory as well. 

Particularly it means that in theories using the principle of local equilibrium 
the entropy function depends only on equilibrium state variables and the entropy 
current is classical. 

In an internal variable theory we use the principle of local state of Kcstin |28U29j . 
implicitly or explicitly supposing that 

- physical quantities defined in non-equilibrium thermodynamics of discrete 
systems and the relations between them are valid in a continuum theory, 
too. 

Particularly it means that in theories using the principle of local state the entropy 
function can depend on internal variables beyond the equilibrium state variables 
and the entropy current is classical. In this way we can get equations incorpo- 
rating higher order time derivatives and we can preserve the stability structure of 
thermodynamics. 

Let us remark that in the literature both concepts appear in a controversial 
treatment (including the mentioned references). For example de Groot and Mazur 
formulates the concept of "local equilibrium" as follows: 

"It will be now assumed that, although the total system is not in equilibrium, 
there exists within small mass elements a state of " local" equilibrium, for which the 
local entropy s is the same function previously defined (for equilibrium systems) of 
u, v and Ck as in real equilibrium." |19| 

The controversy is contained in the concept of equilibrium. A correct formulation 
and understanding incorporates a suitable nonequilibrium theory of discrete (ho- 
mogeneous) thermodynamic bodies including internal variables started by Onsager 
and Machlup [23 EH and continued by others [23 El EI E1EI1 ESI EH , because 
in the above statement an equilibrium system is necessarily time dependent. 

In nonlocal theories we arc beyond the validity of local equilibrium and local 
state and we deal with continuum systems that do not have a discrete counterpart. 
In the following we introduce new kind of internal variables to characterize the 
nonlocality of the continuum. The nonlocal internal variables are current-like in 
some sense, and the extension of the state space appears not only through the 
modified entropy function but also through a modified entropy current. The local 
character of the theory is preserved because the entropy function and the entropy 
current function are local, do not depend on space derivatives. 

The necessity to use a generalized entropy current was pointed out and em- 
phasized several times by Miiller [SHI EH| ■ Regarding the internal variables it was 
Verhas who developed Miillcr's generalization and introduced higher order correc- 
tions with the help of transport of internal variables @U1- The importance of the 
entropy current describing nonlocal effects was pointed out by several previously 
mentioned authors jSHI ETJ EE| ■ 

Our approach will be called weakly nonlocal irreversible thermodynamics because 
we give conditions to transform the entropy inequality into a solvable, quadratic 
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form. In this way the constitutive theory can fully preserve the predictive character 
of irreversible thermodynamics. In this case the constitutive theory is particularly 
easy, one can identify force-current systems with classical, well established meth- 
ods ^3 GI3- ^he applicability of the approach is demonstrated by deriving sev- 
eral particular classical examples like weakly nonlocal heat conduction, generalized 
Ginzburg-Landau and Cahn-Hilliard type equations. 

In the next section we use Liu procedure with a precise prescription of the state 
space, that is free from several weaknesses of the traditional treatment. In this way 
the background and the validity of our approach is clear and we will see among oth- 
ers why and in what sense Ginzburg-Landau-like and Cahn-Hilliard-like equations 
play a distinguished role in weakly nonlocal irreversible thermodynamics. In the 
third section we treat more examples and some other problems by a more heuristic 
approach to show the easy applicability. The last section contains a summary and 
discussion. 

2. Weakly nonlocal irreversible thermodynamics 

There are some basic problems in the heuristic approach of irreversible thermo- 
dynamics that should be addressed in a consistent treatment: 

- One should decide what are the constitutive functions and what are the 
variables (derivatives) they are depending on. That is, we should fix the 
variables and the level of the approximation before the calculations, we 
cannot do it on the fly. 

- A given, fixed state space and constitutive quantities still does not warrant 
that the entropy production can be written in a special quadratic form, as a 
sum of products of constitutive and fixed functions of the state space. One 
should investigate and clarify what are the conditions of that force-current 
structure. We will see, that in classical treatments (classical irreversible 
thermodynamics, classical extended irreversible thermodynamics) the form 
of the entropy current function is an important condition and we will inves- 
tigate whether there are other choices in the considered particular systems. 

In this section we treat the above questions in a systematic treatment We accept 
that the entropy is an extensive quantity, therefore there is an entropy balance and 
we postulate that the Second Law ensures a nonnegative entropy production. Wc 
will determine what are the state variables and what are the constitutive quantities 
in a thermodynamic theory. E.g. the entropy and the entropy current are always 
constitutive, their evolution is determined by the dynamics of other physical quan- 
tities. Our task is to get restrictions that a positive entropy production and the 
other known restrictions (e.g. balance like dynamics) results in the form of consti- 
tutive functions, hence in the dynamics of the physical quantities. We apply the 
Liu theorem, given shortly in a convenient form in Appendix A. First we treat clas- 
sical theories and we generalize them extending the constitutive spaces to include 
nonlocal effects. 

We suppose that the continuum is characterized by physical quantities given by 
variables defined on the basic state space Z and therefore a local basic state of the 
continuum is given by z G Z. The basic state space is called in a slightly different 
context as the space of wanted fields |43j . It is convenient and sufficiently general to 
assume that the basic state space is a finite-dimensional vector space. The processes 
are space-time functions on the basic state space. Therefore in a nonrelativistic 
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treatment, for a given observer a process is z p : R x R 3 — + Z, (i, x) i— > z p (i, x). 
We are to give the dynamic equations determining the processes, therefore we omit 
the treatment of the structure of the process space, because it is determined by the 
dynamic equations. The processes are supposed to be continuously diffcrcntiable 
as many times as it is necessary. Moreover, we are fully aware of the fact that 
the structure of space-time is more refined also in the nonrclativistic case [441 145j , 
but for the sake of simplicity and easier understanding we restrict ourselves to the 
level of treatment usual in the literature of the field. One should be careful in 
applying our following considerations to mechanical theories, because mechanics is 
strongly connected to the structure of space-time |46| . To emphasize formally the 
necessary care, the balances are written in a local frame, relatively to an inertial 
observer. Without any further considerations they are valid in a continuum at rest 
respectively to that observer. 

Wc postulate that the entropy is an extensive quantity with non negative pro- 
duction. This can be formulated in the form of a local balance equation 

(1) d t s + V-3s = <t. >0. 

dt denotes the partial derivative with respect to time. The entropy density s and 
the entropy current j s are considered as constitutive quantities, hence the nonncga- 
tive entropy production o~ s have to be a consequence of pure material properties. We 
assume that the above inequality expresses a pure material property and is fulfilled 
independently on the initial and boundary conditions of the governing equations 
(see [2S1 ) ■ We look for constitutive relations with the help of Liu's theorem to 
ensure that nonnegativity . 

In addition of the entropy and the entropy current there can be other consti- 
tutive quantities in the theory. The domain of the constitutive functions is called 
constitutive state space. The particular theories are determined basically by the 
choice of the basic state space and the constitutive state space. 

2.1. Classical irreversible thermodynamics. In classical irreversible thermo- 
dynamics the basic state space Zqit is spanned by the specific extensive variables. 
Therefore a local basic state of the continuum is given as an element of this space 
a = (a±, aa, a n ) € Zcit, where n is the number of considered extensive quanti- 
ties. The balance expressing the extensivity of the variables is written as follows 

(2) d t a + V • j Q = a a . 

Here the current density of the extensives j a = (ji, J2, j«) is a direct prod- 
uct vector of the conductive current densities respectively and is a constitutive 
quantity. Every current has a tensorial order one degree higher than the corre- 
sponding specific extensive quantity. The dot between the nabla and the current 
denotes a contraction for that extra tensorial order. The constitutive space of clas- 
sical irreversible thermodynamics is spanned by the extensives and their gradients 
Ccit = Zcit x Lin(ZciT^ 3 )- A constitutive state of the material is an element 
of this space (a, Va) e Ccit- The so called process direction space is spanned by 
the next space and time derivatives of the basic state: (dta., <9(Va, V 2 a), where V 2 
denotes the second derivative with respect to the space variables. The constitutive 
quantities that depend on the space-time through the constitutive space are the en- 
tropy, the entropy current and the current density of the extensives (s, j s , j a ). The 
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subset of the constitutive space where they are defined is the constitutive domain 
D C Ccit of the classical irreversible thermodynamic material. 

After these preliminary nominations we are looking for restrictions of the consti- 
tutive functions posed by the balance form dynamic equation and the nonnegativity 
of the entropy' production. Considering the previously given constitutive depen- 
dencies we can write |j2J and as 

(3) d t a + dxia • Va + <9 2 j a : V 2 a - a a = 0, 

(4) d lS d t a + 9 2 s9 t Va + <9ij s • Va + <9 2 j s : V 2 a > 0. 

Here d\ denotes the partial derivative of the constitutive quantities with respect 
to a, by the first member of the constitutive space. <9 2 is the derivative with respect 
to the second member, by Va. Now we are looking for constitutive functions that 
will ensure the validity of the entropy inequality for all possible process directions 
enabled by the balances. 

ApplYING Liu's theorem to the equations above (see Appendix A) the process 
direction space spanned by (d t a, d t Va, V 2 a) corresponds to V. The space into 
which the balances are mapping, that is essentially the basic state space Zqit 
divided by the one dimensional vector space of vectorized time |44| . corresponds to 
V. The application of the theorem results in that there exist T S Zq IT and the 
related Liu-equations are 

(5) dis-T = 0, 

(6) 8 2 s = 0, 

(7) (9 2 j s - rd 2 j a ) s = 0. 

The superscript s denotes the symmetric part of the corresponding second order 
tensor. Moreover, the dissipation inequality will be 

d!j s ■ Va - Tdija • Va + To- a > 0. 

From the first two Liu-equations one can see immediately that the entropy is a 
function of the basic state a and the Lagrange multiplier can be identified with the 
intensive variable defined by Ds = d\s = T. A solution of J7J) results in that 

(8) j s (a,Va) =r(a)j a (a,Va)+j (a), 

where jo is an arbitrary function. From these consequences of the Liu equations 
one can get the dissipation inequality as 

(9) j a -vr + v-jo>o. 

If o~ a — and jo = then the entropy current and the entropy production 
reduces to the well known classical form of irreversible thermodynamics 

js rj a 

and 

j Q (a, Va) • Vr > 0. 

The general solution of the above inequality gives the classical Onsagerian con- 
ductivity relations (see Appendix C), where one can identify the current densities 
of extensives as thermodynamic currents and the gradients of the intensives as 
thermodynamic forces of classical irreversible thermodynamics. 
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with positive semidefinite conductivity matrix L. In general the conductivity matrix 
can depend on the constitutive state and is called as quasilinear. In case of constant 
conductivity one can say that the theory is strictly linear |47j . 

It is interesting to observe that the validity of classical force-current relations is 
connected to the condition that in the entropy inequality a a T should be nonnegativc 
separately to the other parts of the inequality. That condition is not true in general 
and is used to separate dissipative and nondissipative parts of the source term 
<T (For example in a thermomechanics of point-masses one can distinct between 
damping, conservative and exiting momentum sources, force fields.) To avoid the 
details of that problem we demanded that a a = 0. The systems corresponding to 
the above special case are sometimes called as pure dissipative and play a certain 
role in connection to variational principles of irreversible thermodynamics [481 149j . 

Furthermore let us observe that the general entropy current and the source differs 
from the classical expressions in the term jo- Several authors suggest an additive 
supplement to the classical entropy current (e.g. the celebrated K vector of Miiller 
|38| ) . In the following in connection to more general theories we will see that 
such term can have a physical significance which also cannot excluded here (e.g. 
thermodynamic theories of chemical reactions). 

We have seen that in our calculations the classical form of the entropy production 
(the first term in the inequality above) and the classical form of the entropy current 
were consequences of the Second Law, without referring to any 'local equilibrium' 
arguments. Formulating the possible functional dependencies of the corresponding 
quantities and the usual (sometimes) implicit postulates (balance structure, exis- 
tence of entropy and the entropy inequality), we encountered a clear distinction 
between postulates and consequences and the limits of validity of the local equilib- 
rium approximation. For example one can see, that the second half of the principle 
of local equilibrium is completely unnecessary, the validity of the equilibrium state 
functions and the classical form of the entropy current are both consequences of 
structural requirements. 

The choice of the constitutive space looks like a crucial point in the above deriva- 
tion. Remarkable, that an extension of the constitutive space with higher order 
space derivatives does not alter the previous scheme. A short calculation shows 
that in the final Onsagerian conductivity equations the conductivity coefficients 
can be functions of the full constitutive space, that is they can depend on all the 
higher order space derivatives, but without further assumptions the form of the 
entropy production does not change. What can make a difference is the extension 
of the basic state space as we will see in the next subsection. 

2.2. Classical extended irreversible thermodynamics. Extended thermody- 
namics is a nonequilibrium thermodynamic theory, where the currents of traditional 
extensive quantities are introduced as internal variables. The theory is sometimes 
called extended irreversible thermodynamics [501 I51j . extended rational thermody- 
namics [23 or wave approach to thermodynamics [ZJ depending on the details of the 
approach and the school of thcrmodynamicians. These details are not too important 
in our treatment. One of the advantages of using currents as independent variables 
is that in this case the final equations can be compared directly with results of the 
kinetic theory, with the moment series expansion of the Boltzmann equation. This 
comparison gives a remarkable insight into the material characteristics of gases and 
(partially) of fluids. 



WEAKLY NONLOCAL IRREVERSIBLE THERMODYNAMICS 



9 



In the following we comply the procedure of the previous subsection in a some- 
what shortened form. In our treatment the difference is that the basic state space 
of classical irreversible thermodynamics is supplemented by the space of current 
density of the extensives, it is spanned by the variables (a, j a ) £ Zet- Therefore, 
an appropriate constitutive space will be Cet — Z x Lin(Z, R 3 ) x Lin(Z, M 3 ) and a 
constitutive state of the extended irreversible material is (a, j a , Va) £ Cet- In this 
way we supposed that there should be an evolution equation for the current and we 
are to find a possible form of this equation, according to thermodynamic require- 
ments. The conditions are the balances of extensives and the nonnegative entropy 
production. Therefore the balance J5J and the inequality are supplemented by 
the next general form of the searched evolution equation 

(10) d t ja + g = o. 

The constitutive quantities defined on the constitutive space are s,j s and Q. 
Considering the constitutive dependencies, the balances and the entropy inequality 
appear in the form 

d t a + V • j a = cr Q , 

dtja + G = 0, 

d 1 sd t a. + d 2 s ■ d t j a + d 3 s ■ d t Va + 9J S • Va + <9 2 j s : Vj Q + d 3 j s - : V 2 a > 0. 

Here, according to our concise notation, the different dots arc denoting the con- 
tractions of the spacelike tensorial orders (not involved in Zet)- Therefore, because 
of Liu's theorem there exist Lagrange multipliers F e and F„ of the first two equa- 
tions so that Liu equations are 

d lS -r e = o, 
d 2 s-r n = o, 
d 3 s = o, 
9 2 j s -r e i = o, 

(d 3 is) s = 0. 
The dissipation inequality in our case is 

(11) aj. • Va - T n g - T e cj a > 0. 

A solution of Liu equations gives that s and j s do not depend on the gradient of a 
and the partial derivatives of the entropy give the Lagrange multipliers, _Ds(a, j a ) = 
(dis,d2s) = (L e ,L„). Thus the Liu equations give a condition for the entropy 
current as 

(12) d 2 j s (a, j ) = d 1 s(a,j a )I, 

where / is the unit tensor of the basic state space. The postulate of pure dissipa- 
tivity (cr a = 0) also simplify the dissipation inequality and we get that 

(13) 9ij s - Va-d 2 s-£>0. 

However, we cannot give a general solution of (|12|) and 1)13(1 without any fur- 
ther ado. Fortunately we can go farther, because according to theories of ex- 
tended thermodynamics one can exploit some physical requirements writing our 
non-equilibrium entropy function in the following form 

(14) s(a,j ) = s (a) - -j a • m(a,j a ) • j a , 
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where the matrix m is called the matrix of the thermodynamic inductivities, so (a) = 
s(a, 0) is the equilibrium entropy, that depends only on the extensives. It is worth 
remarking that this form, with state dependent thermodynamic inductivities is not 
an approximation if we accept the following requirements: 

- The current variables j a are nonequilibrium, dynamic variables (dynamic 
degrees of freedom) that identically vanish in (local) thermodynamic equi- 
librium. There are extensives with such currents (heat and diffusion cur- 
rents) but there are also evident counterexamples (pressure) . The treatment 
of nonlocality connected to non dynamic variables is beyond the scope of 
this paper. 

- The entropy is a concave function. According to that requirement it is 
concave also in the nonequilibrium part of the state space therefore the 
symmetric part of m is assumed to be positive definite. 

- The derivative of the entropy function characterizes the equilibrium. 

With this premises the mean value theorem results in (| 1 4|) as the most general 
form of the entropy (sec Appendix B). In case of general internal variables where the 
physical meaning is fixed only after the thermodynamic restrictions are applied, the 
Morse lemma can give a more restricted form of the entropy function, with constant 
thermodynamic inductivities |7l I52| . However, in extended thermodynamics the 
current densities cannot be transformed arbitrarily, therefore a constant inductivity 
matrix would be an approximation. 

The form (|14(l of non-equilibrium entropy was introduced first by Machlup and 
Onsager for discrete systems using the time derivatives of the original a quantities 
as dynamic variables It was Gyarmati who suggested the above form for 

field theories using the current densities as dynamic variables |47] . This kind of 
entropy function was originally a characteristic constituent of the wave approach 
of thermodynamics, the other theories realized the principle of local state through 
a generalization of the Gibbs relation. 

The previously introduced non-equilibrium intensives are the following 

nz\ p / • \ _ ds _ ds 1 . dm . 

(15) l e (a, Ja j :- ^--^ 2 Ja 'aa"' Ja! 

ds ( Dm \ 

(16) r„(a,j a ) := — = - I m s + -j a • — J • j Q =: • j Q , 

where the superscript s denotes the symmetric part of the corresponding quantity 
and m is a nonequilibrium inductivity. m is symmetric if m is constant. 

The form (|14f> is still too general to solve l|12|) and l|13fl without further assump- 
tions, therefore we restrict ourselves to the usual approximation, when the matrix 
of thermodynamic inductivities is constant: 

s(a, j Q ) = s (a) - -j • m • j„. 

In this case the solution of (|12fl is similar to that we have got in the previous 
subsection 

j s (aj ) = <9isj a + j (a). 
The dissipation inequality becomes 



{g -m s + \7d lS )-j a >0. 
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We have got a product of fixed functions of the basic state and constitutive 
quantities to be determined. The general solution of this kind of inequalities is 
treated shortly in Appendix C. Here we have got constitutive expressions for Q as 
Onsagerian conductivity relations. Therefore, the form of the evolution equation 
of the current densities is determined by 

j„ =L(g-m s + V9is). 

Here L is a positive semidefinite conductivity matrix, depending on the consti- 
tutive state. Substituting into I|1U|) we get 

dtj a + m s - 1 (L-%-Vd 1 s) = 0. 

2.3. Weakly nonlocal extended irreversible thermodynamics. Our next ex- 
ample is the first weakly nonlocal extension of extended irreversible thermodynam- 
ics. The basic state space is the same as in the last subsection but we will consider 
an extension of the constitutive space. Here we need the first and second deriva- 
tives of the current densities, too. Hence the constitutive space of nonlocal extended 
thermodynamics is C NE t = Z x Lin(Z, M 3 ) x Lin(Z, M 3 ) x Lin(Z, R 6 ) x Lin(Z, WL 9 ) 
and the constitutive functions will depend on the quantities (a, j a , Va, Vj a , V 2 j a ) £ 
Cnet- The balance of the extensives the evolution equation of the currents 
(f 1 Oft and the entropy inequality are to be solved. 

In classical irreversible thermodynamics and extended irreversible thermodynam- 
ics the choice of the state space made possible to prove the local equilibrium and 
local state hypothesis: the entropy function did not depend on the nonlocal part of 
the constitutive space, on the gradient of the extensive quantities. In our case, the 
above extension of the state space does not enable such kind of proof, therefore we 
will postulate it. This is an assumption, we will call that the hypothesis of relocal- 
izability. However, instead of directly prescribing that the entropy does not depend 
on the gradients we make a seemingly weaker hypothesis requiring that the gradi- 
ent of the current evolution equation (|1L)|I does not appear among the constraints. 
Furthermore, we treat pure dissipative systems with dynamic current densities. 

The constitutive quantities defined on the constitutive space are the same as 
previously s, j s and Q. Accordingly, the balances and the entropy inequality appear 
in the form 

d t a + V -ia = 0, 
dtia + G = 0, 

d 1 sd t a + d 2 s ■ d t j a + d 3 s ■ d t Va + <9 4 s : <9 t Vj Q + d 5 s- : <9 t V 2 j Q + 
<9ij s • Va + <9 2 j s : Vj Q + <9 3 j s : V 2 a + <9 4 j s - : V 2 j« + <9 5 j s " V 3 j a > 0. 

Introducing the Lagrange multipliers T e and T n for the first two equalities, the 
Liu equations are 

d lS -T e = 0, 
d 2 s-T n = 0, 

dis = 0, i = 3, 4 

(d 5S y = o, 

(dijsY = 0, i = 3,5. 
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Here the superscript s of dsj s can be expressed best with indexes that (95)yfe(j s )i 
is symmetric in ijk. Because of the larger constitutive space the dissipation inequal- 
ity is considerably longer 

(17) cy s • Va + {d 2 j s - TJ) : Vj a + d 4 j s - : V 2 j a - T n Q > 0. 

A solution of the Liu equations gives that the entropy is a function of a and j a 
solely and does not depend on the gradients at all (that is s(a, j a )) and j s does not 
depend on the gradient of a and the second gradient of j a (that is j s (a, j Q , Vj a ))- 
The partial derivatives of the entropy give the Lagrange multipliers respectively. 
The dissipation inequality reduces to 

v-j s -r e vj„ -r n g>o. 

We can see that the solution of the Liu equations does not give a soluble form 
of the dissipation inequality. However, it is easy to transform it to the traditional 
force-current form after some further, not too restrictive premises. Let us turn our 
attention to the entropy current. A natural physical assumption is that entropy, 
being connected to other extensives, does not flow if the extensive quantities do 
not. Here we should aply the same assumption as in case of classical extended 
thermodynamics suppose the the current of the extensives is a dynamic variable 
and vanishes in local thermodynamic equilibrium. Therefore we can assume that 
the entropy current is zero if the currents of the extensive quantities arc identically 
zero j s (a, 0, 0) = 0. Considering this condition and using the mean value theorem, 
one can get the next functional form 

(18) j s (a,j a , Vj a ) = B(a,j a , Vj a ) • j a , 

where the current intensity factor B, can be supposed to be continuously diffcrcn- 
tiable in a neighborhood of the local equilibrium state. We should remark again 
that this form itself is not an approximation, it is general with the above conditions. 
One can observe that the tensorial order of the current intensity factor B is one 
order higher than the corresponding current. As we mentioned above, the classical 
form of the entropy current was investigated and generalized by several authors. 
The above kind of miltiplicative generalization is due to Nyfri [^Sj ■ 
Now the entropy production can be written as 

(19) (B-r e /):Vj a + (V-B + £m)-j a >0. 

where m is the inductivity function defined in l|16(l . The general solution of the 
inequality is 

V • B - dtja ■ m = Luj a + Li 2 Vj Q , 

(20) B-TJ = L 21 j a + L 22 Vj a . 
We can eliminate the current intensity factor and get 

(21) m* • dtja - V ■ (L 3X j fl + L 22 Vj„ + TJ) = - (L n j a + L 12 Vj a ) . 

The star denotes the transpose. One can see, that contrary to the common 
postulate |39j in extended thermodynamics, this evolution equation does not have 
a balance form in general. It could be transformed formally into a balance form 
with special source in case of constant inductivitics. 

Let us remark that a smaller extension of the constitutive space with only the 
gradient of the current densities would not have been enough to get a gradient 
dependent entropy current. 
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2.4. Weakly nonlocal internal variables. In the previous subsections the phys- 
ical meaning of the internal variables was fixed, they were the current densities of 
the thermodynamic extensives. However, it is a restriction, the usage and applica- 
bility of internal variables is far more broader |54l I55j . We can meet them under 
different names as order parameters, dynamic degrees of freedom, dynamic variables 
in different physical theories with a slightly different meaning. They arc applied 
characterizing materials with complex microstructure, properties of gases near crit- 
ical states, etc. The extent and method of applying thermodynamic restrictions in 
internal variable theories can be very different in different approaches. In thermo- 
dynamic theories the application of internal variables is connected to the principle 
of local state and therefore to discrete systems [21] ■ Their dynamics is restricted 
or generated constitutively by the Second Law and results in ordinary differential 
equations. In discrete systems, i.e. in homogeneous bodies, the formalism is well 
known and frequently applied, however sometimes the dynamic aspects are hidden 
(a good summary of the dynamics of the simplest discrete thermodynamic system 
without internal variables is given from several point of view in |57U58ll^H| and with 
internal variables in |33U35j ). On the other hand, as we mentioned in the introduc- 
tion, these internal variable theories can be considered as a "relocalization" of some 
memory, inertial effects of the equilibrium variables. The homogeneous dynamics 
(ordinary differential equations) in continuous systems has some important disad- 
vantages: the local state seems to be too tight from several points of view. First 
of all some physical phenomena cannot be modelled with them (e.g. thickness of 
shear bands, structure of surfaces, ...). Moreover, the numerical solution of coupled 
ordinary and partial differential equations depends on the discretization, on the 
mesh. 

The situation is entirely different in case of weakly nonlocal equations. Here one 
can distinguish between essentially two independent approaches. One of them is 
based on the fundamental balance equations introducing gradient dependent terms 
directly into the constitutive functions, which are generated using the local equilib- 
rium or the local state hypothesis (see e.g. |59j and numerous examples in rheology 
and mechanics in general |60M13lH4*] l. 

The other approach is even more intuitive. There one assumes that the local 
entropy function depends on the gradients of the corresponding physical quantities 
(usually on a rather simple way) and a kind of 'nonlocal intensives' are generated 
by the functional derivatives of the entropy function. This method is used to get 
the Ginzburg-Landau equation for a so called 'non conservative' physical quantity 
(internal variable), and the Cahn-Hilliard equation for a 'conservative', extensive 
one (model A and B with the terminology of Hohcnberg and Halperin [HI])- The 
general equation for the noncquilibrium reversible- irreversible coupling (GENERIC) 
is based on a generalization of this kind of approach |62l . 

The relation of the two approaches is not clear at all. With the traditional 
approach (based on local equilibrium and state) one cannot get the most important 
nonlocal theories of contemporary physics (e.g. Ginzburg-Landau) on the other 
hand the variational derivations are too intuitive and they have nothing to do with 
the balances of the fundamental physical quantities. Criticizing the mentioned 
derivations of Ginzburg-Landau and Cahn-Hilliard equations, Gurtin emphasizes 
the importance of the separation of balances from the constitutive properties: 
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"My view is that while derivations of the form ... are useful and important, they 
should not be regarded as basic, rather as precursors of more complete theories. 
While variational derivations often point the way toward a correct statement of 
basic laws, to me such derivations obscure the fundamental nature of balance laws 
in any general framework that includes dissipation. " |64| 

Gurtin himself uses a method exploiting his " microforce balance principle" , 
which is a new principle again, and we can apply Gurtin's further critical remarks 
used against the variational approach also to his propositions: 

- The derivations limit the manner how other terms (e.g. rate terms) enter 
the equations. 

- The derivations are specific, it is not clear how they can be generalized for 
general thermodynamic processes including not only mechanics, diffusion 
or the dynamics of the order parameters near critical states. 

In this section we will see how far can we go with relocalized theories, with 
the methods of irreversible thermodynamics. We introduce nonlocal internal vari- 
ables connected to the entropy current. First we will see the weakly nonlocal 
irreversible thermodynamic theory of pure dynamic internal variables without fur- 
ther constraints, after that we investigate the case of a nonlocal extension of a 
thermodynamic theory containing also extensive variables. 

2.4.1. Weakly nonlocal dynamic internal variable - Ginzburg- Landau equation. Our 
task is to find weakly nonlocal evolution equation of an internal variable that cor- 
responds to the requirement of nonnegative entropy production. The basic state 
space, Zi, is spanned by the vector denoting an array of internal variables with 
different tensorial character. We are to find an evolution equation of the internal 
variable in the form: 

(22) d t £ + T = 0. 

The constitutive space is spanned by £ and its first and second gradients: (£, V£, 
V 2 £). Therefore C7 = Z^xIm(^,M 3 )xLm(^,R 6 ) is the constitutive space of the 
nonlocal dynamics of an arbitrary internal variable. The constitutive quantities are 
the entropy, the entropy current and the form of the evolution equation (s,j s , J 7 ). 
Therefore, the positive entropy production is supplemented by 11221) in the Liu 
procedure. 

dt£ + T = 0, 

dxsdti + d 2 s • 9 t V| + 8 3 s : d t \7 2 £ + 0J 8 • V£ + d 2 j s : V 2 £ + 8 3 j s - : V 3 <£ > 0. 

Here we applied the hypothesis of relocalizability, as in case of extended ther- 
modynamic systems. According to Liu's theorem there exist a T, to be determined 
from the Liu equations, which can be written in a particularly simple form 

d lS T = 0, 

d 2 s = 0, 

d 3 s = 0, 

d 3 j s = 0. 



The dissipation inequality in our case is 
(23) c\j s • V£ + 9 2 j s : V 2 ^ - TT > 0. 



WEAKLY NONLOCAL IRREVERSIBLE THERMODYNAMICS 



15 



The solution of Liu equations gives that s depends only on the internal variable 
£, r = Ds(£) is the derivative of the entropy and j s does not depend on the 
second gradient of Unfortunately these considerations do not simplify the entropy 
inequality at all, wc should look for additional conditions. As previously we assume 
that ^ is a dynamic variables. Therefore there is no flow of entropy and the derivative 
of the entropy is zero if the internal variables are zero, that is j s (0,0) = and 
r(0) = 0. The first condition is quite natural. The second condition requires that 
the zero derivatives characterize the equilibrium of the material. Thus we can write 
the entropy current in the following form 

j s (£,V£)=A(£,V£)r(£). 

This form of the entropy current is general with the assumptions above. As 
previously, we assume, that A is a continuously diffcrcntiablc function. In this case 
the entropy inequality turns out to have the form 

(v ■ a - .F)r + a • vr > 0. 

Let us remark, that in (|23(l there are two constitutive quantities (the entropy current 
and J-) and three additive terms. To simplify the inequality for constructing a force 
current system there is no other choice that we have done here: we should unite 
two of the terms with some reasonable physical assumption. Now we can solve the 
resulted inequality and get an Onsagerian structure. 

(24) V-A-.F = L n r + L 12 vr, 

(25) a = L 21 r + L 22 vr. 

Here L is the positive semidefinite constitutive conductivity matrix. Eliminating 
A from (JUJ) and we get 

(26) d t £ + V • (L 21 r + L 22 Vr) = L n L + L 12 VL. 

This is a balance equation with a special current and source term. To recognize 
the structure of the equation let us consider a single scalar internal variable £ and 
isotropic material. In this case there is no cross effect between the thermodynamic 
forces of different tensorial order Ln = and L 2 i = 0, and the diagonal elements 
of the conductivity tensor L are scalars L\\ = l\ and L 22 = Z 2 / (according to 
Curie principle or more properly the representation theorems of isotropic tensors). 
Therefore equation 1)26(1 simplifies to 

(27) d£ = hT - V ■ (J a Vr c ). 

Here we have got an equation that is similar to the well known Ginzburg-Landau 
equation. However, there are some differences. 

In the second term of the right hand side under the space derivatives is T 
instead of £, contrary to the Ginzburg-Landau equation. The consequence 
of this difference is that the equation has the traditional homogeneous ther- 
modynamic equilibrium solutions characterized by the condition L = 0. 

- The sign of the material coefficients is determined by the Second Law, not 
by direct stability considerations. 

- It is straightforward to extend the derivation to nonlinear and anisotropic 
cases. However, the nonlincarities and anisotropics show a different struc- 
ture than in the original Ginzburg-Landau equation. 
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In the following (|27J) will be called thermodynamic Ginzburg- Landau equation 
|f)5j . This equation was derived by Verhas under some slightly different assump- 
tions, as a governing equation for the transport of dynamic degrees of freedom 

mi nam 

Remark 2.1. In all of the previous considerations the only mathematical step that 
seems to be restrictive from a physical point of view is the (continuous) differen- 
tiability of the current intensity factor. That condition can be weakened, but it will 
exclude the appearance of terms like jo in case of classical and extended thermo- 
dynamics. The question is not mathematical and academic at all. We should not 
forget that the question of differentiability plays an important role in constitutive 
modelling and in thermodynamics in general. For example some phase boundaries 
are classified according to the differentiability of equilibrium state functions. More- 
over, recently an experimental indication of phase boundaries was found also in the 
non- equilibrium state space, connected to internal variables. That concept makes 
possible the unification and the elaboration of a common thermodynamic frame for 
theories of failure and for theories of elasticity in mechanical continuum (e.g. in 
case of damaged brittle materials |f)7l I68j ). 

2.4.2. Weakly nonlocal dynamic internal variable and extensive variable - Cahn- 
Hilliard equation. In our last example there is an extensive thermodynamic variable 
and also an internal one that is not the current density of the extensive as in ex- 
tended irreversible thermodynamics. Therefore the basic state space of the physical 
quantities consists of the extensive quantity a with a balance like evolution equa- 
tion J5J and a dynamic internal variable £ with an evolution equation (|22J) . whose 
form is not specified. We are looking for nonlocal extension of the internal variable, 
therefore the constitutive space consists of the constitutive spaces of classical irre- 
versible thermodynamics and the previous pure nonlocal internal variable theory 
together C IE = Z a x x Lin(Z a , M 3 ) x Lm(%M 3 ) x Lm(Z f ,l 6 ). Our consti- 
tutive quantities are the entropy, the entropy current, the current of the extensives 
and the evolution equation of the internal variable (s, j s , j Q , T). All these quantities 
are functions defined on the constitutive space spanned by (a, £, Va, V£, V 2 £). We 
are looking for the restrictions on the constitutive quantities imposed by the non- 
negative entropy production JQ) and the balance J2J with zero source term (pure 
dissipativity) and the condition of relocalizability. Considering the later we do not 
impose the space derivative of the evolution equation of the internal variables as 
an additional constraint. Because of the constitutive dependencies, the constraints 
and the entropy inequality can be written in the following form 



d t a + dij a ■ Va + d 2 ] a ■ V£ + d 3 j a : V 2 a + d 4 j a : V 2 £ + d 5 j a - : V 3 £ = 0, 

d t i + T = 0, 

d lS d t a + d 2 sd t £ + d 3 s ■ d t Va + d 4 s • t V£ + d 5 s : <9 t V 2 £+ 
SJ, • Va + d 2 j s ■ V£ + 6> 3 j s : V 2 a + d 4 j s : V 2 £ + d 5 j s - : V 3 £ > 0. 
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According to Liu's theorem there exist functions L e and to be determined 
from Liu equations written as 

<9is - T e = 0, 

d 2 s - r 6 = o, 

dis = 0, i = 3,4, 
(d 5 s) s = 0, 

(d 3 j s - r e 9 3 j Q ) s = o, 
(d 5 j s - r e 9 5 j a ) s = o. 

Here the superscript s of the last equation can be expressed best with indexes that 
(c?5)r/fc(j)z is symmetric in ijk. A solution of the Liu equations results in that the 
entropy depends only on the extensive and on the internal variable. The derivative 
of the entropy gives the Lagrange multipliers Ds(a, £) = (d a s,d^s) = (L e ,T^). A 
consequence of the last two equations is 

(28) j 8 = r e j +jo(a,£,VO- 

Here jo is an arbitrary function of the variables denoted above. In the first term 
we can recognize the classical entropy current. This solution of the Liu equations 
simplifies the dissipation inequality to 

(29) Vr e -j a + V-j O -L e ^>0. 

This form is simple but unfortunately we cannot solve it withoutany forther ado. 
We require that in case of zero internal variable the entropy flow reduces to the 
classical form. The dynamic property of the internal variable gives the additive 
term of the entropy current in the following form 

jo(a,£,V£) = A(a,£,V£)IY 

Now the dissipation inequality transforms to 

j a • VL e + A • VL 5 + (V • A - T)^ > 0. 

As in the previous sections wc reckon the entropy as a primary quantity. There- 
fore we arrived to a force-current system, to a sum of products of undetermined 
and given functions of the constitutive state. Therefore the general solution the 
inequality results in the Onsagerian conductivity equations in a quasilinear form, 
with a positive semidefinite conductivity matrix L, depending on the constitutive 
state 

(30) j a = Lnvr + Li 2 vr, + l 13 t ( , 

(31) a = L2ivr + L 22 vr c + L 23 r<r, 

(32) V-A-T = LnVL + Li 2 Vr e + L 3 r e . 

Substituting (JSOJl into © and (j3H> into O we arrive the following system of 
transport like equations with a definite balance form 

(33) a t a+ v- (Lnvr + L 12 vr s + L 13 r e ) = 0, 

(34) 04 + v- (L 21 vr + L 22 VL 5 + L 23 r e ) = LnVr + L 12 vr f + l 3 iy 

On the other hand this choice of the basic state space indicates that the above 
set of evolution equations can be considered as generalization of the Cahn-Hilliard 
equation. To see the connection one must eliminate the internal variable. Let 
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us consider a very special example, where our basic state space is spanned by a 
single scalar extensive a and a single scalar internal variable £. Furthermore, let be 
the material isotropic and all the material coefficients are constant (strictly linear 
approximation). Finally, let us write the entropy in the canonical form (|14f) . but 
with constant inductivities m. Now there arc two coupled terms in the linear laws, 
because the current of the extensives j a and the current multiplier A have same 
tensorial order 

(35) j Q = Z u Vr - Z 12 mV£, 

(36) A = l 21 WT - Z 22 ?™V£, 

(37) V-A-J" = -l 3 m£. 

Here the matrix of the conductivity coefficients is positive definite to ensure 
nonncgativc entropy production (h > 0, In > 0, I22 > and Zn/22 — ^12^21 > 0). 
We can eliminate V£ from the first two equations to get an expression of the current 
multiplier A as: 

A = Wialja + (hi - l22luhi)VT. 
We can put this formula into (|37|) and so eliminating £ from (j35(l we get 

(38) j Q = luV(r + aAT) + /3VV • j Q + h^dtt 

where a = I11I3 (I12I21 — ^22^11) and (3 = l^ 1 h2- we can neglect the last term in 
the above expression if I12 = or the change of £ is very slow. So we eliminate the 
internal variable and got (0) and (|38(l as a system of equations to be solved. After 
some straightforward manipulation we finally get 

(39) d t a + lnA(T -aAT) - (3Ad t a = 0. 

In our case, when a represents a single scalar physical quantity (e.g. mass den- 
sity) then In, ol and /? are positive numbers. Equation l|39|l is similar to the 
Cahn-Hilliard equation, but not the same, therefore it will be called thermodynamic 
Cahn-Hilliard equation. It is worth to give a short comparison. 

- The last extra term does not appear in the traditional equation. Gurtin 
derived a similar additional term with the help of microforce balance (a 
strong additional assumption) in the special case when the variable a is 
the mass density |64j . As regards some generalizations and the properties 
of the solution see (SHI E0I- This term is similar that one could get in 
case of the Guyer-Krumhansl equation and there are several experimental 
indications for this kind of generalization with the Cahn-Hilliard equation, 
too |71j . Supposing that Z22 = 0, this additional term vanishes. Let us 
recognize that this rate dependent term appears in the linear Onsagerian 
conductivity equations at the same thermodynamic approximation level as 
all the others. 

- Instead of the extensive a, the intensive variable T is under the second 
Laplace operator, contrary to the original Cahn-Hilliard equation. The 
equilibrium solutions of the thermodynamic Cahn-Hilliard equation are ho- 
mogeneous. 

- Every material parameter in (|39|l is positive, as a consequence of the Second 
Law, a nonnegative entropy production. In the original equation the signs 
of the coefficients are fixed according to stability considerations seemingly 
independently of the entropy balance. 
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- The nonlinearities and anisotropies show a different structure than in the 
original equation. 

In most of the previous derivations of the Cahn-Hilliard equation, a " gener- 
alized chemical potential" appeared by analogy, and was inserted without 
any further ado into the mass balance. Here in our derivation we have got 
a real generalization of the diffusion current in the same sense and with the 
same method as in the original irreversible thermodynamic approach. 

- Treating the derivation of the thermodynamic Cahn-Hilliard equation a 
scalar extensive and a scalar internal variable was introduced. As a can 
denote the Cartesian product of several extensives, the coupling in conduc- 
tivity equations (|35() - (|37|1 can become more involved, even in the simplest 
isotropic case. On the other hand, the introduction of an internal variable 
with the tensorial order of the current (one order higher than the extensive) 
results in the very same structure of the final Cahn-Hilliard type equation. 

3. Heuristic weakly nonlocal irreversible thermodynamics 

In the previous section weakly nonlocal thermodynamics was treated in a more 
or less mathematical level. Our point of view, the leading idea behind the pre- 
vious section was to give the background and the conditions of applicability of 
the classical irreversible thermodynamic approach far beyond the validity of local 
equilibrium. We have seen, that the essence of the classical method of irreversible 
thermodynamics, the construction of proper force-current systems can be applied 
not only to systems with inertia, but for essentially nonlocal systems, too. More- 
over, the local character of the theory can be preserved introducing special internal 
variables, current multipliers. 

The general basic conditions to create solvable force-current systems were the 
following: 

(1) Entropy - there exist an extensive entropy with nonnegative production. 
Entropy and entropy current are constitutive functions. The entropy is the 
primary constitutive function of any thermodynamic theory. 

(2) Pure dissipative systems - the source term is zero in the balances of consid- 
ered extensives. 

(3) Relocalizability - the entropy, as fundamental potential function expressing 
the properties of ordinary (equilibrium thermodynamic) systems does not 
depend on gradients of the basic variables. 

In weakly nonlocal theories this is a postulate, in case of more restricted 
constitutive spaces is a consequence of the thermodynamic requirements as 
we have seen in the previous sections. 

(4) Dynamic internal variables - are those that vanish in (local) equilibrium. 
This assumption results in two consequences for sufficiently smooth consti- 
tutive functions: 

If the entropy is a concave function of the internal variables (in the 
noncquilibrium part of the state space) then we arrive at the Gyarmati 
form of noncquilibrium entropy (|14|> - 
- The Nyfri form generalized entropy current 1)18(1 . 
Here, in this section, based on the above properties, we treat some further weakly 
nonlocal relocalizable systems in a more heuristic way. This heuristic point of view 
shows, how one can get solvable (Liu) equations and solvable entropy inequality 
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in the more exact and formal mathematical procedure. In the first subsection we 
investigate the parade ground of thermodynamic theories and develop weakly non- 
local heat conduction. After that we show that a current multiplier of current 
multipliers can have an importance and derive Gurtin's term in Ginzburg-Landau 
equation. Finally we recall that current multipliers arc not unique, there are sev- 
eral possibilities to transform the dissipation inequality into a solvable form. Here 
we investigate extended thermodynamics as a special internal variable theory and 
investigate the conditions to get a balance like evolution equation for the internal 
variables. 

3.1. Weakly nonlocal heat conduction - Guyer-Krumhansl equation |72l 

173) . The different special theories of thermodynamic origin as diffusion, mechanical, 
electromagnetic interactions, etc. all can show phenomena that can be explained 
by weakly nonlocal extensions of the classical equations. Let us investigate more 
closely the most frequently treated basic example of all thermodynamic investiga- 
tions, the phenomena of heat conduction. Now the single extensive quantity is the 
internal energy, therefore a = (ai) = u. Instead of citing the previously derived 
formulas we will apply a simplified, more heuristic procedure. The example demon- 
strates the role of the continuization hypotheses of local equilibrium, local state and 
rclocalizability giving the different heat conduction theories respectively. We will 
see that the Guyer-Krumhansl equation is the thermodynamically consistent first 
weakly nonlocal extension of the Cattaneo-Vernotte equation and the structure of 
the equations is independent of the particular material properties. 

The balance of internal energy is the starting point. Let us denote the heat 
current density by j a = (j u =)q and we assume that the system is purely dissipative. 



(40) d t u + V • q = 0. 

For the case of simplicity, we consider isotropic material and constant thermal 
inductivity m. In this case m = ml, therefore the entropy function l|14|) can be 
written as 

(41) s(ji,q) = .s (u)-imq 2 , 

where m > 0. The Nyi'ri form of the entropy current is 

(42) j s = B„ ■ q. 

Now the entropy production (|19f) can be calculated by composite derivations 

dts + V • j s = 
d u sd t u + d q sd t q + V • (B„ • q) = 

(43) (B, - -I) : Vq + (V • B„ - md t q) • q > 0. 

We can recognize in our special case. According to Curie principle (rep- 
resentation theorems of isotropic tensors, e.g. |74l 175) ) there is no cross coupling 
between the different thermodynamic interactions. Moreover in 1)20(1 Ln = II, 
where I is a scalar and L22 can be written with indices as 

{L22)ijki = hSikSji + hSuSjk + hSijSki- 
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Therefore the linear approximation of Onsager H2()|) reduces to 



q = Z(V • B M - md t q) 




This is exactly the Guyer-Krumhansl equation for the heat current, introduced 
to describe the thermal properties of some crystals at low temperatures ■ Orig- 
inally it was derived using kinetic physics and later two- fluid hydrodynamics [7flj . 
At the last decades there were several attempts to get the equation from pure 
non-equilibrium thermodynamics |77l 1781 1731 181)1 15T1 1411 152] . However, all of these 
derivations contain several ad- hoc assumptions, whose generalization for more dif- 
ficult situations can be very ponderous if not impossible (e.g. two-fluid hydrody- 
namics is an internal variable theory from a thermodynamic point of view, but with 
a particular interpretation of the internal variables) . One of these assumptions is a 
balance like dynamics for the heat current q that is a result of the general structure 
in our case. If we do not eliminate B„ from the material equations, then we can 
get a dynamic equation of the heat current that has a special balance form (more 
properly it will be an equation of Ginzburg-Landau type) indicating that the cur- 
rent intensity factor can be considered as a current density of the heat current q 
|4(Jj . We analyze further that property later. 

In the special case of h = fo = h = 0, we get the Cattaneo-Vernotte equation for 
the heat current. The conductivity equations indicate clearly that the local state 
hypothesis is applied, the entropy current has its traditional form j s = j g /T, but the 
entropy itself depends on the dynamic noncquilibrium variable q. The Cattaneo- 
Vernotte equation is a wave like equation for the heat conduction phenomena. It 
was Gyarmati who first derived it in the frame of irreversible thermodynamics with 
the help of entropy function l|41[l |4"7] . Further, we get the Fourier heat conduction 
equation when m = 0, in local equilibrium, when the entropy depends only on 
internal energy u. 

In case of constant material coefficients (m, I, l\, ...) it is easy to derive a telegraph 
type equation directly for the temperature T. We should introduce the specific heat 
u = cT and substitute l|40|l and its time derivative into the divergence of (|45|l 



This reduction is remarkable and discussed from several points of view. First 
of all in several practical problems the material parameters cannot be considered 
as constant quantities (see e.g. [HSJ )■ On the other hand the reduction of a vector 
and a scalar equation (|4(J|) and (|45|l into a single scalar one results in a loss of 
information. The system l|4U|) and (|45|l is equivalent to (|46|l if the heat current field 
is rotation free. 

3.2. Nonlocal internal variables. 



(46) 



ImcduT + cd t T + I A- + lc(h + h + h)Ad t T = 0. 



3.2.1. Generalized thermodynamic Ginzburg-Landau equation. In this subsection 
we develop the results of subsection 12.4.11 and show how can an additional rate 
term appear in the Ginzburg-Landau equation. 
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As previously our task is to find an evolution equation of an internal variable 
that corresponds to the requirement of nonnegative entropy production. The basic 
state space is the vector (moreover Banach) space of the variable £ and we do not 
know a dynamic equation for the processes. We suppose a rclocalizable theory, 
therefore the entropy function depends only on the internal variable. We will use 
the notation := Ds(£). If £ = then 1^(0) = 0, because £ is a dynamic variable. 
As regards the entropy current we apply the previous physical assumptions also in 
this case: if £ was zero then there is no entropy flow. Moreover, in the light of 
this assumption and according to the mean value theorem, the entropy current 
can be written as a linear function of the derivative of the entropy, as in classical 
irreversible thermodynamics j s = AF^. Here we introduced a nonlocal internal 
variable, a current multiplier A. We can go farther, and introduce the current 
multiplier into the basic state space. For the sake of simplicity, we suppose that 
the entropy function does not depend on A, that is we are not interested in the 
associated memory effects, we arc investigating only the nonlocal extension. The 
form of the entropy current is similar to that of the Cahn-Hilliard equation: 

j s = Ar e + B- A. 

Here B is a current multiplier of the first current multiplier A. This form ex- 
presses the physical assumption, that there is no entropy flow when the internal 
variables £ and A are zero. The tensorial order of A is one order higher than 
the tensorial order of £ and the tensorial order of B is one order higher than the 
tensorial order of A as it should be considering the most general linear connection. 
Therefore the entropy production follows as 

(47) dtsfc) + v • i s = (dt£ + v • (A)r 5 + a • (vr ? + V • B) + B : VA > 0. 

It is straightforward to put down the Onsagerian conductivity equations, but in 
general one cannot simplify them. Therefore, we will treat here only the simplest 
situation, when £ = £ is scalar, the material is isotropic and the approximation is 
strictly linear (the conductivity coefficients are constants). Now the conductivity 
equations are reduced to the following form 

(48) dtt + V-A = 

(49) A = Z 2 (Vr c + V-B) 

(50) B = ZgVA + Z|(VA)* + ZfV • A, 

where Zi, I2, 13, Z§, Z| are positive, scalar, constant coefficients. Now a simple calcu- 
lation eliminates A and B from the above equations: 

(51) a t £ = hT £ - Z 2 A(1 + z x z 3 )r £ + z 3 Ad t e 

Here I3 = Z3 + Z3 + Z3. The last term, that is additional to (|27() corresponds to the 
generalized Ginzburg-Landau equation of Gurtin. One can get back the thermody- 
namic Ginzburg-Landau equation i(27|) in the special case when = Z| = Z| = 0, 
that is B = 0. The positivity (positive definitcness in a more general situation) 
of the material coefficients is ensured by the Second Law. However, we should ob- 
serve, that this generalization has not changed the characteristic thermodynamic 
term which differs from the original Ginzburg-Landau form: F^ stands under the 
Laplace operator instead of £. 

We can continue the introduction of new nonlocal internal variables, putting 
B into the basic state space. In this case B becomes an internal variable and 
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we can introduce a corresponding current multiplier. Continuing this procedure, 
we can develop a whole phenomenological hierarchy of weakly nonlocal transport 
equations of higher and higher orders. The further research in this direction has a 
special importance for the kinetic theories. The outlined phenomenological hierar- 
chy of nonlocal equations can suggest an approach similar to the momentum series 
expansion |84l I8()| . On the other hand it is straightforward to extend the above 
treatment (with Liu procedure or without it) considering memory effects, too. 

3.2.2. Are the dynamic equations of internal variables balances or not? In section 
12 .31 we have seen that in extended irreversible thermodynamics, when the currents 
of the extensives are the internal variables, the final constitutive-evolution equation 
of the currents did not have a balance form in general. Later, in section 12.4. II for 
a pure internal variable the derived evolution equation l|34[l had a definite balance 
form. The reason of the difference is hidden in the way of the generalization of 
the entropy current. The classical form reads as a product of the currents of the 
extensives and the corresponding entropic intensives 

(52) j s = d a sja 

In the generalization of Nyfri, in section l2~3l devcloping weakly nonlocal extended 
irreversible thermodynamics, the entropy current is assumed to be a homogeneous 
linear function of the currents of the extensives 

(53) j,=B-j a . 

On the other hand, in case of dynamic variables that vanish in equilibrium and 
the entropy characterizes the equilibrium the entropy current was assumed to be a 
linear function of the entropic intensives, conjugated to the corresponding extensive 
variable. 

(54) j s = A ■ d s s. 

This form is a different generalization of (|52|l than i|53|) and was applied in sec- 
tion l2.3l relatcd to model B, leading to the thermodynamic Cahn-Hilliard equation. 
Extended thermodynamics looks like a special case of the system treated in 12.31 
because currents of extensives can be considered as special internal variables. How 
the two assumptions are related to each other? Could we arrive to an evolution 
equation with balance form also in case of extended thermodynamics as it is gener- 
ally expected [SHUSH? To answer this question we should investigate the forms of 
the entropy inequality more closely. Instead of simply refcring to the results of the 
previous sections we will give again a less refined, briefer, more direct, althought 
sloppier, traditional irreversible thermodynamic approach. 

To get the entropy production of weakly nonlocal extended thermodynamics <|19[) 
we consider a system with a basic state spanned by extensive a and internal variable 
£ together. In case of extended thermodynamics the (dynamic) internal variable is 
the (dynamic) current of the extensives £ = j a . Assuming a local entropy function 
s(a, ja) one can get lfH3|) by substituting (JHHJ), ifTH)) and (0) into the entropy balance 
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as follows 

d t s + Vj 6 = 
d a s ■ d t a + djs ■ d t j a + V • (B • j„) = 
~d a s ■ V • j a - djs ■ Q + B : Vj„ + V ■ B • j a = 

(55) (B-d a sI) : Vj Q + (V-B + em)-j a > 0. 

Here we exploited the Gyarmati form entropy l)14|l and wrote djs = — m • j a . 
Because the entropy function is the primary constitutive quantity we can conclude 
that H55fl resulted in a solvable inequality for he two (!) constitutive functions B 
and Q. 

On the other hand, when the basic state space is the same as previously but 
with a general dynamic internal variable one evidently cannot assume any more 
that the entropy current is proportional to the dynamic variable and one cannot 
write the entropy current in a one term form of 1)54(1. However, a splitting into two 
terms, a classical and a non classical one, works well. Therefore we will assume the 
following configuration of the entropy current 

(56) j s = d a s) a + B' 

In this case a we can go forward similarly to the previous derivation and substi- 
tute J5J) and Gyarmati entropy into the balance of entropy to get 

d aS ■ dtEL + d iS ■ d t i + v • {d a s3a + b' • £) = 

~d a s • V • j a - d e s ■ T + d a sV ■ j a + j a ■ V9 a s + B' : V£ + V • B' ■ £ = 

(57) B' : V£ + j a • V<9 a s + (V • B' + mf) •£•>(). 

This a solvable inequality for he three (!) constitutive functions B', T and j a . 
Moreover it reduces to the previous (|55|l one in the special case of £ = j a , when 
Q = T and B = B' + <9 a s/. Therefore j s = (B' + d a sl) j a . However, the solution of 
the inequality does not result in a balance form evolution equation for the internal 
variable, similarly to the extended thermodynamic case. 

In section ?2 . 4 . 21 the entropy current was assumed to have the following, different 
form 

(58) j s = d a sj a + A • d^s. 

Now, we can substitute J2J and l(28|l into the balance of entropy to get the entropy 
inequality in the following configuration 

d t s + V-j s = 

d a s ■ d t 8L + d^s ■ dti + V • {d a sj a + A ■ 9js) = 

~d a s ■ V • j - <%s ■ T + d a sV ■ j a + j a ■ V9 tt s + A : V<%s + V • A • d^s = 

(59) A : V9 e s + j a • Vd a s + (V • A - T) ■ d f s. > 0. 

This a solvable inequality for he three (!) constitutive functions A, T and j Q 
again. Moreover, as we have seen in section 12.4.21 the solution results in a balance 
form evolution equation of the internal variable without applying Gyarmati entropy 

In the special case of extended thermodynamics, when £ = j a the entropy in- 
equality transforms to 

<7 S = A : Vdjs + j a ■ Vd a s + (V • A - T) ■ djs > 0. 
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As one can see the above inequality does not have a general solution being a sum 
of three product terms with two constitutive functions. However, with Gyarmati 
entropy djS = — m ■ j a and the inequality further transforms to 

a s = m • A • Vj a + (V • (8 a al - m ■ A) - m • F) ■ j a > 0. 

Here we can recognize the classical solvable extended thermodynamic form with 
B' = m • A and the entropy current in the Nyfri form as j s = B ■ j a = (— m ■ 
A + d a sl) • j a . The inequality was transformed into a solvable force-current system, 
but the balance form of the arising evolution equation of the internal variable is 
destroyed. 

As a summary of the previous investigations we have seen that in case of dy- 
namic internal variables in general the balance form of the evolution equation is a 
consequence of natural physical conditions but in case of extended thermodynamics 
the limited constitutive structure of the theory seem to destroy the balance form 
of the evolution equations as general solutions of the dissipation inequality. 

Let us remark here, that the easier formal applicability and heuristic power of 
the classical irreversible thermodynamic method can support but cannot replace 
the deeper constitutive investigations of the first section. First of all because the 
postulates of the more exact method were clearer and weaker than the heuristic one 
(e.g. the form of the entropy current (|28J) was consequence of the Liu equations). 
On the other hand one should identify correctly the basic state space and the con- 
stitutive space to see beyond the rather formal treatment of current-force systems 
of irreversible thermodynamics. 

4. Summary, conclusions and discussion 

In the previous sections we have seen that it is possible to extend irreversible 
thermodynamics to nonlocal phenomena with the help of generalized internal vari- 
ables. The extension is based on the generalization of the entropy current. Partial 
differential equations were constructed constitutively as evolution equations of the 
internal variables by a force-current structure. The structure of nonlinearities, the 
sign of coefficients were determined by the requirement of nonnegative entropy pro- 
duction. We have investigated the thermodynamic background of several classical 
equations of mathematical physics containing higher order space derivatives. 

In the first section we have seen that the choice of the variables, basic and 
constitutive state spaces are crucial steps in thermodynamic theory construction. 
Basic state spaces are spanned by the physical quantities for that one wants to 
find an evolution equation. Constitutive state spaces are spanned by variables of 
the basic state space and their derivatives, the ones, the constitutive functions are 
supposed to depend on. The extension of the constitutive state space determines 
whether an evolution equation of an extensive or an internal variable considers 
nonlocal (or memory) effects or not. 

In this paper we investigated the nonlocal extension only of dynamic internal 
variables. According to this view we can call a basic internal state variable nonlocal 
if the constitutive space contains at least its second order space derivatives. The 
following table summarizes the hierarchy of continuum theories from that point of 
view, a denotes the array of extensive variables and £ an array of internal variables 
as above. 
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Basic state 


Constitutive state 






Local 


Nonlocal 


Classical 


a 


(a,Va) 




Extended 


(a,jo) 


(a,j ,Va) 


(a,j a ,Va,Vj a ,V 2 j a ) 


Internal variable 






(£,v^,v 2 £) 


Classical + internal 
variable 


(a,€) 


(a,e,Va) 


(a,£,Va,V|,V 2 £) 



- The first row of the table concerns classical irreversible thermodynamics 
where the basic state space is spanned by specific extensive quantities. 
Classical irreversible thermodynamics is a local theory, in the sense of local 
equilibrium. Investigating its nonlocal extensions one can easily conclude 
that in case of relocalizable theories (when entropy depends only on the 
basic state) , the structure of the equations does not change and only the 
classical constitutive functions become gradient dependent. 

In the second row, where the basic state space is spanned by the exten- 
sive variables a and their (dynamic) current densities j a one get classical 
extended thermodynamics. Extending the constitutive space with the gra- 
dients of the current densities gives a nonlocal theory and leads to Guyer- 
Krumhansl-like equations. 

- For pure weakly nonlocal dynamic internal variables £ the nonlocal ex- 
tension of the classical ordinary differential equations results in Ginzburg- 
Landau-like equations. 

- An extensive variable with a weakly nonlocal dynamic internal variable 
leads to Cahn-Hilliard-like equations. 

Our investigations were restricted to relocalizable theories, where the nonlocality 
can be grasped through nonclassical internal variables, introduced into the entropy 
current. The examples can be continued, other kind of basic and constitutive states 
and couplings can lead to other weakly nonlocal thermodynamic compatible equa- 
tions. However, we may not forget that there are genuine weakly nonlocal theories 
and equations, with large importance and ample experimental evidence. A typical 
example is the traditional, non thermodynamic Ginzburg-Landau equation. For- 
tunately our methods are applicable to embed also true weakly nonlocal equations 
into noncquilibrium thermodynamics. The key aspects of Liu's procedure regard- 
ing that question and the first investigations related to the compatibility of the 
Ginzburg-Landau equation with noncquilibrium thermodynamics are given in |85j . 

We can observe that the emerging classification scheme is more refined that was 
given by Penrose and Five (EH|, Hohenberg and Halperin j^J or Gurtin j^l]. The 
investigation of several existing weakly nonlocal equations (e.g. different phase field 
equations, the complex Ginzburg-Landau or Kardar-Parisi-Zhang) from the point 
of view of nonequilibrium thermodynamics can reveal richer structures, new terms, 
clearer conditions, surprising interrelations, hence a broader range of applicability 
than usually considered. 

In this paper it was shown that contrary to the common belief, irreversible ther- 
modynamics can be extended far beyond the usual local equilibrium or local state 
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hypotheses, one can consider dissipation related to nonlocal effects, too. A rigorous 
treatment showed that the seemingly intuitive steps of irreversible thermodynamic 
modelling arc well supported and explained with Liu procedure. Several interesting 
consequences, new terms in the traditional equations and possible sources of phys- 
ical generalizations were explored. There arc several experimental and theoretical 
evidences that the characteristic additional term, containing the time derivative of 
the Laplacian of the basic state variable (A<9 t ), appears in all of the treated classi- 
cal equations (Guyer-Krumhansl, Ginzburg-Landau and Cahn-Hilliard). That term 
cannot be derived from variational considerations without any further ado, because 
its Frechet (strong) derivative is not symmetric |87| . 

An other important remark, that the stability structure of nonequilibrium ther- 
modynamics is preserved in the weakly nonlocal extension. Therefore the entropy 
can serve as an Liapunov functional with appropriate initial and boundary condi- 
tions. These possible boundary conditions arc more general than that were given by 
Penrose and Fife in their investigations regarding the thermodynamic compatibility 
of Ginzburg-Landau and Cahn-Hilliard equations. The physical difference is that 
now Liapunov functions can be generated for several different open systems. This 
kind of generalization goes farther than it was possible in homogeneous systems 
(see e.g. |37p and also requires a specification of the function spaces where the 
solutions of the differential equations are investigated. 
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6. Appendix A: Liu-procedure 

The Liu-procedure is based on an application of a theorem in linear algebra. That 
theorem can be considered as a consequence of the Farkas' lemma [SS], well known in 
linear programming. More properly it is a special case of the so called affine Farkas' 
lemma [EH] (that can be easily derived from the original Farkas' lemma [HD])- The 
recognition of this connection between the Liu theorem and Farkas lemma is due 
to some recent researches of Kirchncr and Hauser |91j . 

In this appendix we formulate and prove that theorem of Liu |22j in a simple 
and general form best suited for our purpose. The proof below is essentially that 
of Liu. 

V and V are finite dimensional vector spaces, V* and V* denote their duals. 

Theorem 6.1. (Liu) A e Lin(V, V), B e V, a e V* and b e R. TTiera a-x+6 > 
holds for si! x £ V such that A • x + B = Oy */ and only if there is a A 6 V'* such 
that 

(60) a- A- A = Oy, 

(61) 6-A-B > 0. 

f6'0|) is called Liu-equation and \61)) is the dissipation inequality in the thermo- 
dynamic literature. 



Proof: 



28 



P. VAN 



The backward direction is easier. If there is a A with the above property, then 
for all x G V, < (a - A • A) • x - A • B + b = (a • x + b) - A ■ (A • x + B) and the 
statement follows. 

Now let us prove the forward direction. We introduce the notation Aff = 
{x e V | A ■ x + B = Oy } ■ First of all we should see that KerA C Kera. Really 
if y G KerA and x G Aff then for all r £l,x + ry£ Aff, therefore according to 
the conditions of the theorem < a ■ (x + ry) + b = a ■ x + b + ra ■ y. For arbitrary 
r the inequality is true if and only if a • y = therefore KerA C Kera. 

For the annullators of the kernels 

(Kera) 1 - C (KerA) 1 -, 

because Vx G KerA C Kera, therefore, if p G (Kera) 1 - then p ■ x = 0. 

On the other hand KerA* = (RanA)- 1 , because if Vx e V y • A ■ x = = 
x • A* ■ y then from the left hand side of the equality y G (RanA) 1 - and from the 
right hand side y G KerA*. For finite dimensional V, A** = A and (A- 1 )- 1 = A, 
therefore we can transform the equality into the form 

(KerA) 1 - = RanA*. 

Putting together all these statements we can see, that a G (Kera) 1 - C (KerA) 1 - = 
RanA*, therefore there is a A G V* such that a = A* - A. In finite dimension follows 
(|6(J[I . Therefore we can write for all x with A • x + B = Oy that < a • x + b = 
A- A-x + & = 6-A-B. 
The proof is complete. 

Remark 6.1. If A is surjective (dim(RankA) = dimY' ) then A is unique. 

7. Appendix B: Mean- value theorem 

We applied Lagrange's mean value theorem several times in the paper. Here we 
formulate it in a way that is sufficiently general for our purpose. One can find a 
proof for example in |92|. 

Theorem 7.1. (Lagrange) Let X be a normed space. F : X — > R, x, y G X. If F 

is differ entiable on the closed section [x, y], then there is a z G [x, y], so that 

F(y)-F(x)=DF(z)-(y-x), 
where DF(z) denotes the derivative of F at z. 

8. Appendix C: Solution of the Dissipation Inequality 

In a force-current form the dissipation inequality can be written as a sum of the 
product of constitutive quantities and variables in the constitutive space. In this 
case a solution can be given in the form of the (nonlinear) Onsagerian conductivity 
relations, we formulate the statement and give the simple proof following Gurtin 

EH- 

Theorem 8.1. Let V be a finite dimensional vector space. J : V — > V*. x i— > J(x) 
is a continuously differentiable function. If J(x) • x > for all x G V then J(x) = 
L(x)x where L : V -> Bilin(Y). 

Proof: 

From the inequality follows that for all positive number A and x G V, J(Ax)-Ax > 
0, therefore J(Ax) ■ x > 0. In the limit when A goes to zero we get that J(0) ■ x > 0. 
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That can be true for all x G U if and only if J(0) = 0. The statement follows 
according to the mean value theorem. 

Remark 8.1. The structure is not unique in case of nonlinear relations. 
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